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At present it remains to address why the fast solar wind is fast and the slow wind is slow. Recently we have shown that the field line 
curvature may substantially influence the wind speed v, thereby offering an explanation for the Arge et al. finding that v depends on 
more than just the flow tube expansion factor. Here we show by extensive numerical examples that the correlation between v and field 
line curvature is valid for rather general base boundary conditions and for rather general heating functions. Furthermore, the eff'ect of 
field line curvature is even more pronounced when the proton-alpha particle speed difference is examined. We suggest that any solar 
wind model has to take into account the field line shape for any quantitative analysis to be made. 
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Half a century after its prediction IT] and direct detection ||2l, 
how the solar wind is accelerated to its highly supersonic ter- 
minal speed remains elusive. Among the extensive exper- 
imental findings that provide stringent constraints on these 
mechanisms, those concerning the in situ wind speed v are 
of primary concern from both theoretical and practical stand- 
points: V and the magnetic field polarity near the Earth are 
the two main parameters involved for predicting geomagnetic 
activities, a major player in space weather ID |4] |5] 01 ■ A 
widely used technique to predict v is the Wang-Sheeley-Arge 
model |7|, which is based on two well-established statisti- 
cal relations between v and the properties of the flow tube X 
along which the wind flows. First, v is inversely correlated 
with the coronal expansion rate /c of X IS). Second, for a 
given /c, V is inversely correlated with the angular distance Ot, 
of the footpoint of X to its nearest coronal hole boundary |f7l|. 
Given the importance of the wind speed prediction, the two 
results deserve a sound physical explanation. While the for- 
mer is well supported by a number of wave/turbulence-based 
solar wind models (see Q and references therein), the latter 
has so far attracted Uttle attention. 



Close inspections of coronal images indicate that the fur- 
ther away the field line footpoint is from the coronal hole 
boundary, the more curved the field line is (see e.g., the com- 
posite image Fig.l in ifTOl ). Consequently, the effects of Oi, 
on the wind speed may be explained in terms of field line 
curvature. A detailed modeling study using a sophisticated 
Alfvenic-turbulence-based heating mechanism suggests that 
it is indeed so ( ifTTI . hereafter paper I). The model distin- 
guishes between electrons and protons, adopts reasonably 
complete energy equations including radiative losses, elec- 
tron heat conduction, and solar wind heating. It also employs 
the Radiation Equilibrium Boundary (REB) at the base which 
is set at the mid-Transition region, whereby the base electron 
density automatically adjusts to the downward heat flux from 
the corona 11211131 . Compared with a straight flow tube with 
identical radial distribution of the magnetic field strength, a 
curved tube produces a significantly reduced flow speed at 
1 AU (by up to ~ 30%), accompanied by an enhanced mass 
flux density. This effect was interpreted, in view of the gen- 
eral framework by lfT4l . as a result of the enhanced energy 
deposition in the near-Sun region in the curved case relative 
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to the straight one since the elementary reduction in the wave 
action flux density is a factor of dl/dr larger (Eq.(7) in paper 
I), where / is the arclength and r the heliocentric distance. 

The aim of the present paper is to extend paper I by ex- 
amining whether the effect of field line curvature discov- 
ered there is robust enough: Is it vaUd only for this particu- 
lar boundary condition or this particular heating mechanism? 
Will the inclusion of alpha particles, the second most abun- 
dant ion species, invalidate the conclusion? The latter is of 
particular importance since alpha particles, with their mass 
being 4 times the proton one, may play an important role in 
regulating the wind parameters ifTSl IT6l fTTl . This paper is or- 
ganized as follows. We will give a brief overview of the solar 
wind model in section[T] then present the numerical results in 
section[2] Finally, section [3] summarizes the results. 

1 Model Description 

Let us start with considering a time-independent solar wind 
that consists of electrons (e), protons (p), and alpha parti- 
cles (a). Each species s {s = e, p, a) is characterized by its 
mass nis, electric charge e.,, number density n.,, mass density 
Ps - rtsnis, velocity Vj, temperature Ts and partial pressure 
Ps - rtskBTs- Here kg is the Boltzmann constant. We may 
express in units of the electron charge e, i.e., = Z^e 
with Ze s -1 by definition. Assuming axisymmetry, quasi- 
neutrality and quasi-zero current, neglecting solar rotation, 
one may readily decompose the vector equations in the stan- 
dard 5-moment treatment into a force balance equation across 
the meridional magnetic field B and a set of transport equa- 
tions along it (for a formal development, see section 2.1 
in IfTSl ). This transport equation set reads 
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where the subscript s refers to all species {s - e,p,a), while 
k stands for ion species only (k = p, a). The gravitational 
constant is denoted by G, Mq is the mass of the Sun. The 
prime ' denotes the derivative with respect to the arclength 
I, measured from the footpoint of a magnetic field line. The 
momentum and energy exchange rates due to the Coulomb 
collisions of species s with the remaining ones are denoted 
by dMJdt and dE Jdt, respectively. Their expressions may 
be found in (191 for which we use a Coulomb logarithm 
In A = 23. As for a^, it represents the acceleration exerted on 
ion species k. Moreover, qs is the heat flux carried by species 



s, and Qs stands for the heating rate applied to species s from 
non-thermal processes. L is the radiative loss function in- 
volved in the electron energy equation only. It is assumed to 
be optically thin and follows the analytic form given in ll20l . 
To simplify our treatment, we neglect ion heat fluxes (qk - 0), 
and adopt the Spitzer law for the electron one qe - -kT^^^T'^ 
where = 7.8 x 10"^ erg K^^^^ cm"' s"' ||2T1 . 

As detailed in paper I, rather than solving for a fully 
self-consistent multi-dimensional solution, we solve Eqs.([T]i 
to (O on a prescribed magnetic field instead, the reason be- 
ing that in a multi-dimensional model the eff'ects of mag- 
netic field expansion and field line curvature cannot be dis- 
tinguished. In practice, we employ the analytic model 
given by Il22l for which the reference values [Q,K,ai] = 
[1.5, 1, 1.538] are used. The resulting magnetic field con- 
figuration agrees favorably with the coronal images obtained 
with SOHO/LASCO ll23l . and has been widely seen as being 
representative of a minimum solar corona. In line with the 
Ulysses measurements, the magnetic field strength at 1 AU is 
essentially latitude-independent and is 3.5 y Il24ll25l . All the 
presented results will be computed along the field line that 
reaches 1 AU at 88° colatitude, nearly maximizing the eff'ect 
of field line curvature. 

Specifying the boundary condition (BC) at the base r = 
1 Rq is crucial for constructing wind solutions. We will 
test three BCs: REB, TR, and CB. By REB, we mean the 
Radiation Equilibrium Boundary where a base temperature 
Tso = 5 X 10^ K is set, the electron density «eo = 4.8 x lO-'^eo, 
in which the subscript represents the base value lfT2l [T3l . 
Besides, TR (CB) refers to the case where we set the base 
at the lower Transition-Region (coronal base), Tso = 10^ K 
(lO*" K), the corresponding «eo being 10'' (1.5 x 10'^) cm"\ 

The source terms ak and Q, remain to be specified. In our 
study they depend on whether alpha particles are included, 
and on how the base BC is implemented. Electron heating 
is applied for the REB and TR boundaries to compensate for 
the significant radiative losses at the corresponding temper- 
atures. It is in the form = Geoexp(-///£/), with geO (in 
10"^ erg cm^"* g ') and /</ (in Rq) being 1 .7 and 0.06 for REB, 
2.9 and 0.06 for TR. For ion heating, either the wave heat- 
ing scenario or some empirical heating function is employed. 
The former choice is used only in the absence of alpha parti- 
cles. In this case a wave transport equation (Eq.(4) in paper 
I) is included where the parallel propagating Alfven waves 
lose their energy via both doing work on the flow throught 
the ponderomotive force (ak, see Eq.(2) in paper I) and dis- 
sipation at the Kolmogorov rate. (For more recent develop- 
ments on the role of Magnetohydrodynamic waves in ener- 
gizing plasmas, see e.g., 1261 .) We use a base wave ampli- 
tude of 27, 18, and 29 km s"', a base value of correlation 
length being 1.5, 2, and 2x10"* km, for the REB, TR and CB 
boundaries, respectively. These parameters are within obser- 
vational limits, and produce solar wind solutions whose ter- 
minal speed and proton flux agree with in situ measurements. 
When ions are heated in an empirical fashion, we assume the 
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Table 1: Combinations of base boundary conditions, electron and proton heating, and inclusion of alpha particles 
section Base Boundary Condition e heating ion heating a-particles included? 
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See text for details. 
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external energy comes from an ad hoc energy flux dissipated 
at a constant length /ad, the resulting total heating rate be- 
ing Q — FEBI{BEhd)tx^{-llli,i). Here Fe is the input flux 
scaled to the Earth orbit = 1 AU, and Be is the meridional 
magnetic field strength at r^. Q is apportioned between pro- 
tons and alpha particles, and some energy that alpha particles 
receive turns out necessary to be in the form of momentum 
deposition, 

G, ^ Qa Pa 

Qp Pp 

Qa + PaVaCla - Qa, — - , 

Qa PaVl) 



where vo - ^jkBTplnip is some characteristic speed. The 
momentum deposition to protons is neglected. This way for 
empirically heating the ion fluids seems rather involved, but 
actually it manages to capture some essential features of wave 
heating in a multifluid treatment (for details, see ll27l . sec- 
tion 3.1). In practice, we use Fe - 1.4 erg cm"^ s"' and 
'ad = 1 -75 Rq, regardless of the presence of alpha particles, 
such that the resulting proton flux and speed at 1 AU agree 
with the in situ measurements. When alpha particles exist, 
we use Xe - 6 and xd = 2, the purpose being to generate a 
wind profile that corresponds to a positive proton-alpha speed 
difference (see section |23] for further details.) 

Table 1 summarizes our computations. The first col- 
umn indicates in which section the combinations of BCs and 
ion/electron heating in the same row are tested. For bound- 
ary conditions imposed at the base, REB (Radiation Equilib- 
rium Boundary), TR (Transition Region), and CB (Coronal 
Base) conditions are experimented with. For electron heat- 
ing, Y (N) refers to the case where the heating is switched on 
(off). As for the ion heating, "Wave" and "empirical" repre- 
sent the case where the ions are heated by turbulent Alfven 
waves and in an empirical fashion, respectively. The last col- 
umn indicates whether or not alpha particles are included in 
the computations, with Y and N meaning that they are and 
are not, respectively. 




Case 



C REB 
S REB 
C TR 
S TR 
C CB 
S CB 



10 



100 



Figure 1 Effects of field line curvature on wind speed under the influence 
of diff'erent base boundary conditions. Here the proton speed Vp is shown as 
a function of heliocentric distance r. The solid (dashed) curves are for the 
solutions that incorporate (neglect) the field fine curvature. Black, red, and 
green curves represent solutions using the REB, TR, and CB boundaries, re- 
spectively. Moreover, the vertical bars represent wind speeds measured by 
tracking small density inhomogeneities (blobs) in SOHO/LASCO images as 
given by 1281 . 

2 Numerical Results 

Now we are in a position to examine how robust the effects of 
field line curvature are by varying our choices of base BC, the 
proton heating, and the inclusion of alpha particles. For any 
combination of the three choices, we compute two solar wind 
solutions that differ only in whether or not the field line cur- 
vature is considered. If not (yes), the solution is labeled S (C), 
short for straight (curved). Before proceeding, let us briefly 
discuss what we mean by "curvature". Similar to paper I, 
here it is defined in an intuitive rather than a precise manner, 
and refers to how significantly the field line differs from a ra- 
dial shape. Locally speaking, at any given point along a field 
line, the curvature can be precisely defined, say, in terms of 
the curvature radius H. However, the effect of field Une shape 
on the wind parameters turns out to be an integrated one to 
which the shape at each and every point contributes. A non- 
local definition of the field line shape by, say, integrating 1 fR 
is for sure more precise, but evaluating its correlation with the 
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computed wind parameters is beyond the scope of the present 
paper 

2.1 Influences of the base boundary condition (without 
alpha particles) 

Let us for now neglect alpha particles but focus on the in- 
fluence of the boundary condition imposed at the base. Fig- 
ure 1 presents the distribution of the proton speed Vp with the 
heliocentric distance r. It compares the solutions obtained 
in the curved cases (solid lines) with those in the straight 
cases (dashed lines). The black, red, and green curves are for 
REB, TR, and CB, respectively. Moreover, the vertical bars 
represent wind speeds measured by tracking small density 
inhomogeneities (blobs) in SOHO/LASCO images as given 
by ||281 . One can see that all lines, both solid and dashed, 
possess a dip below 2 Rq, indicating that this local minimum 
in Vp results from the drastic, lateral expansion of the flow 
tube rather than its curvature. On the other hand, all the solid 
curves reproduce the SOHO/LASCO measurements equally 
well, and they reach similar values at 1 AU. To be specific, 
the REB, TR, and CB solutions yield a Vp of 333, 355, and 
351 km s respectively. The corresponding values in so- 
lutions with identical BCs but neglecting field line curvature 
are 461, 462, and 482 km s"' respectively. Moreover, the pro- 
ton flux at 1 AU (in 10** cm- s ') is found to be 3.65(3.14) in 
the curved (straight) case for REB, 3.86 (3.47) for TR, 3.57 
(3) for CB. In other words, one finds a reduction in Vp, rel- 
ative to the straight case, of 27.8% for REB, 23.2% for TR, 
and 27.4% for CB, and an enhancement in proton flux by 
16.1% for REB, 11.4% for TR, and 18.9% for CB, respec- 
tively. Hence, the significant reduction in the terminal proton 
speed and enhancement in the proton flux due to the inclu- 
sion of field line curvature is not specific to the particular 
REB condition employed in paper I, but happens for rather 
general boundary conditions. 

2.2 Influences of proton heating 

Is the effect of field line curvature present only for the par- 
ticular ion heating mechanism used in paper I? This is exam- 
ined in Figure 2 where we again neglect alpha particles but do 
employ a different ion heating function, namely the empirical 
heating rates detailed in section[T] Figure 2 presents the radial 
distributions of the proton speed Vp, for both cases C (solid 
lines) and S (dashed lines). Comparison of Fig. 2 with Fig.l 
indicates that the proton speed profile in Case C is steeper 
than in the wave-driven model, nonetheless the effect of field 
line curvature remains the same. To be specific, the termi- 
nal speed is reduced from 583 in the straight to 399 km s 'in 
the curved case. Likewise, the proton flux (in 10** cm^ s ') 
increases from 2.61 to 3.47. This change in the proton param- 
eters can be readily explained using the same argument as in 
the wave-driven case. The energy deposition in the subsonic 
region can be expressed as an energy flux density scaled to 



1 AU, Fsub = L QBe/ B(dl/dl)dr, evaluating which results 
in ^£[1 - exp(-/c//ad)]- Here the subscript C represents the 
values at the critical point where Vp equals ^ksiTf. + Tp)/mp. 
It follows that if in Cases C and S the critical point lies at the 
same radial distance, Fjub is larger in Case C given that Ic 
is larger. And from the general framework presented in lfT4ll 
it can be readily deduced that Case C will yield an enhanced 
proton flux and a reduced proton terminal speed. Actually 
this effect is compounded by the fact that in Case C the criti- 
cal point is located further away from the Sun: rc =4.13 and 
3.51 Rq in Cases C and S, respectively. 



2.3 Influences of alpha particles 




Figure 2 Effects of field line curvature on the proton speed for an empir- 
ically heated wind. The solid (dashed) curves are for the solutions that in- 
corporate (neglect) the field line curvature. The vertical bars have the same 
meaning as in Figure 1. As for the asterisks, they represent the critical points 
where the proton speed equals the sound speed (see text for details). 

Finally let us address the role of alpha particles. For sim- 
plicity we employ the empirical ion heating rates, since we 
have akeady shown that adopting some more sophisticated 
mechanism like wave heating will not change our result as 
far as the effect of field line curvature is concerned. Figure 3 
examines how the flow parameters at 1 AU change with vary- 
ing iiapfi, the alpha abundance at the coronal base. These 
include the proton speed Vp,£, proton flux density {npVp)E, 
and proton-alpha speed difference Vap,E - (va - Vp)E- The 
solid and dashed curves are for Cases C and S, respectively. 
Figure 3a indicates that for riapfl in the range from 3 x 10"^ 
to 0.1, the field line curvature consistently reduces the termi- 
nal proton speed and enhances the proton flux density. For 
^ap.o S 10"^, introducing alpha particles has Uttle effect on 
the proton parameters and the reduction in £ by the field 
line curvatures is ~ 183 km s"', or by ~ 31.5%. On the other 
hand, for nap,o = 0.1, Case C yields a Vp^E of 331 km s"', or 
down by 29.6% relative to Case S where Vp ^ is 470 km s"'. 
Hence whether or not the alpha particles are included, intro- 
ducing field line curvature has the same effect on the proton 
parameters in interplanetary space. 
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Figure 3 Effects of field line curvature on flow parameters for an electron- 
proton-alpha wind. Displayed as a function of the base alpha abundance 
"ff/j.o = (na/np)o are values at 1 AU of (a) the proton speed VpE and pro- 
ton flux density (npVp)E, and (b) the proton-alpha speed difference Vap_£ = 
(Vtt - Vp)£. The solid (dashed) curves are for the solutions that incorporate 
(neglect) the field line curvature. 

The proton-alpha speed difference v^p has been a subject 
of extensive study (e.g., ||29ll30l and references therein). Fig- 
ure 3b shows that the field line curvature significantly reduces 
x'ap- This is true throughout the range of napfi explored. For 
instance, the amount of reduction in Vap,E is 133 km s"' for 
nap.o = 3 X 10"-'* (from 264.3 in Case S to 131.4 km s"' in 
Case C), and is 195.5 km s"' for iiap.o =0.1 (from 200 in 
Case S to 4.54 km s"' in Case C). This significant reduction 
in Vop at large distances turns out to be reminiscent of what 
happens in the corona: although not shown, the alpha tem- 
perature is substantially higher in Case S than in Case C from 
1 .5 Rq onward, resulting in a substantially larger alpha pres- 
sure gradient force and hence a significantly larger proton- 
alpha speed difference. From this we may conclude that, the 
effect of field line curvature should be taken into account for 
a quantitative investigation into the low-latitude slow wind, 
and its role in regulating the ion speed difference is even more 
pronounced. Once this effect is introduced, one may expect 
to see some significant revision to some recent studies (e.g., 
iflTllSTl ). where the flow tube is taken to be radially directed. 

3 Conclusions 

The present study has been motivated by the fact that for 
space weather applications, the near-Earth solar wind speed 
V is of crucial importance and is being routinely predicted 
by the Wang-Sheeley-Arge (WSA) model Q. This model 
is based on the relations between v and the coronal expan- 



sion rate of the solar wind flow tube f8'| as well as the angu- 
lar distance of the tube footpoint to its nearest coronal hole 
boundary |7J • While the former correlation has been exten- 
sively examined, the latter has been given little physical in- 
terpretation, and was only recently explained in terms of the 
field line curvature [1 1 1 1 : the more curved the flow tube (or 
equivalently the magnetic field line) is, the higher the proton 
flux density, and the lower the terminal proton speed. Given 
that the conclusion was established only via a limited num- 
ber of computations, and given that providing a sound phys- 
ical explanation for the WSA model puts considerably more 
confidence in this model, we show in this paper by extensive 
numerical examples that the correlation between the terminal 
proton speed and field line curvature is valid for rather gen- 
eral base boundary conditions and for rather general heating 
functions commonly adopted in the literature. Furthermore, 
its validity is not violated by the introduction of alpha parti- 
cles, the second most abundant ion species in the corona and 
solar wind. An immediate consequence of our new compu- 
tations is that any solar wind model, for the fast and slow 
wind alike, has to take into account the field line shape for 
any quantitative analysis to be made. Beside, the effect of 
field line curvature is even more pronounced when examin- 
ing the proton-alpha particle speed difference, thereby calling 
for the urgent inclusion of this effect in multi-fluid models of 
the low-latitude solar winds where the ion speed difference is 
among the primary concerns. 
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